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A possible supersymmetric solution to the discrepancy between B — > <pK s and 

B — > 7]' Kg CP asymmetries 
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We present a possible supersymmetric solution to the discrepancy between the observed mixing CP 
asymmetries in B — > 4>Ks and B — + rj ' Ks- We show that due to the different parity in the final 
states of these processes, their supersymmetric contributions from the R-sector have an opposite 
sign, which naturally leads to S$k s 7^ Sri 1 Kg' We a ^ so cons ider the proposed mechanisms to solve 
the puzzle of the observed large branching ratio of B — ► rj' K and their impact on S v 'k s . 



Various measurements of CP violation in B factory exper- 
iments have been opening a new era for the phenomenol- 
ogy of supersymmetric (SUSY) models. While in the 
standard model (SM), all the CP violating phenomena 
have to be explained by a single phase in the Cabbibo- 
Kobayashi-Maskawa matrix, the SUSY models include 
additional new sources of CP violation. Since these new 
effects can manifest themselves in the CP asymmetries of 
various B-meson decays, the recent observed large dis- 
crepancy among the CP asymmetries of B — > J/tpKs, 
B — > (pKg and B — > rj'Kg have raised high expectations 
for indirectly unveiling low energy SUSY [1]. 

The measurement of the angle of the unitarity triangle 
(3(4>x) by the so-called golden mode B — > J/ipKg [2,3]: 



S 



J/-4>K S 



sin 2/3(20i ) = 0.734 ±0.054 



(1) 



is in a good agreement with the other measurements 
based on the standard model analysis. Flowingly, it has 
been shown that the effect from the SUSY particles in 
the box diagram which leads to the B° — B mixing is 
typically small [4]. On the contrary, in summer 2002, the 
B factory experiments reported a surprising result for the 
measurement of /3 by using the B — > <j>Ks process. Since 
in the SM, B — > J/ipKs and B — > <j>Ks have the same 
B°-W mixing part and do not have any additional CP 
violating phase in the decay process, the same value of 
sin 2/3 was expected to be extracted from them. Thus, 
the discovered large discrepancy [5,3] 



S. 



<t>K s 



-0.39 ±0.41 



(2) 



has created quite a stir. Several efforts to explain this 
experimental data, in particular, by using SUSY mod- 
els, have been made. In Ref. [6], it has been shown that 
this phenomena can be understood without contradict- 
ing the smallness of the SUSY effect on B — ► J/ipKs 
in the framework of the mass insertion approximation 
which allows us to perform a model independent analy- 
sis of the SUSY breakings [7]. In this approximation, 
SUSY contributions are proportional to the mass in- 
sertions (Sf^AB where i,j and A, B are the generation 
and chirality indices, respectively. While the measure- 
ment of B — > Jji^Ks implies the smallness of (5f 3 )AB, 
(A,B = L,R), the different generation mass insertion 
contributing to the B — ► 4>Ks process, (#23) ab, can devi- 



ate S#k s from S 



J/i>K s 



In this letter, we discuss another 



measurement of sin 2/3 [3,8] 



S. 



■n'Ks 



0.33 ±0.34 



(3) 



which has been thought to be problematic [1]. Since 
B — > ri'Ks gets contributions from (8% 3 )ab, S V 'k s an d 
S^k s were expected to display similar discrepancy from 
Sj/ipK s - We will first show that although the magnitude 
of the SUSY contributions to these processes are indeed 
similar, B — ► rj' Ks has an opposite sign in the coefficient 
for the RL and RR mass insertions, which can naturally 
explain the experimental data. In fact, there is another 
open question on the B — > rj' K process, the observed 
unexpectedly large branching ratio [9]. We will further 
investigate the proposed new mechanisms to enhance the 
branching ratio of B — > rj K and their impacts on S^Ks ■ 
The Effective Hamiltonian for the AB = 1 processes 
induced by gluino exchanges can be expressed as 



H. 



AB = \ 
eff 



Gj_ 

V2 



v tb v t 



2—3—6,g 



ao t + dOi 



(4) 



where the operators Oi can be obtained from Oi by ex- 
changing L <-> R. The Wilson coefficients Ct and Ci are 
proportional to Sll,lr and Srr^rl, respectively. The 
definition of the operators and Wilson coefficients (and 
the effective Wilson coefficients below) can be found in 
[6]. Employing the naive factorisation approximation 
[10], where all the colour factor N is assumed to be 3, 
the amplitude for the B — > cf>K process can be expressed 



A((f>K)-- 



Gf 

'7f 



6 
i=3 



a 



eff 



cf 



(0i?°|O 4 |B°) (5) 



where the matrix element is given by 

= l^x, x, l -x 



(i 



6) (6) 



with X = 2F 1 B ^ K (ml)Um 4> (p K ■ e ). F?^ K (m%) - 
0.35 GeV is the B — K transition form factor and /a = 
0.233 GeV is the decay constant of the <f> meson. Since 
both Fi^ K (m^) and f<f, are insensitive to the chirality 
of the quarks, we used 



(^lO^B ) = (<j>K° 1 0^°) 



1 



to derive Eq. (5). On the other hand, the amplitude for 
B — > rfK can be written by: 



Gi 



v tb v* 



.i=l 

i=3 



V2 

where we used 

(r } 'K \O i \B ) = -{n'K°\d i \B ) 



(v'K^OilB ). (8) 



(9) 



which is derived by the fact that the decay constant of 
rj is sensitive to the chirality of the quarks. The matrix 
element is given by: 



(rjK Q \O x \&) 
(rjK°\0 3 \S°) 
(rj'K°\0 4 \B°) 
(rj'K°\Or,\B°) 
(rj'K°\O 6 \B ) 
with 



-X 2 , (r,'K°\O 2 \B ) = X 2 , 
lx 1 + 2X 2 + ^X 3 , 



X\ + 



:*3, 



?±X 1 -2X 2 + (-l + ^)X 3 , 
R 1 X 1 -^X 2 + (~+R 2 )X 3 



(10) 

(11) 
(12) 

(13) 
(14) 



XI = -{m% m%,)F^{m\ Y-^=f K , 



X v > 

7f' 



X2 = -(>»% *K)/ir^f> 



X3 = -(m| - m 2 K0 )F^ K {ml)^2f%- f^Y v ,, 



Rl 



2m 2 K0 



{mb - md)(m s + md) 



R2 



2(2m 2 KO - ml) 
(mb - m s )(m s + m s ) 



where F^ 3 ~* 1T (q 2 ) = 0.3 GeV is the B — it transition form 
factor and fx^) = 0.16(0.13) GeV is the decay constant 
of K(w) meson, q is the momentum transfer of the b — > s 
transition. X v i = 0.57 and Y v i = 0.82, which correspond 



to 9 n 



-20°, represent the rate of the uu + dd and ss 



component in the rj [11,12]. We use the following quark 
masses, (m dl m s ,m b ) = (0.0076,0.122,4.88) GeV. The 
tree contributions to B — > rj' Ks is found to be less than 
1% and can be ignored. 

Numerical results on the ratio between SM and SUSY 
amplitudes for mg ~ m^ = 500 GeV are obtained as [6] 



<t>K s 



V 'K S 



(0.23 + 0.04z)[(<5i L ) 23 + (4 K ) 23 ] 

+(95 + 14i)[(^ B ) 2 3 + (<5^)23] (15) 
(0.23 + 0.04i) [(8 d LL ) 23 -(SUhs] 

+ (99 + 15z)[(0 23 -(4l)23]. (16) 



where a parameter q 2 is chosen to be m 2 /A. The varia- 
tion of q 2 within the range of m 2 /6 < q 2 < ml/3 causes 



±30% of theoretical uncertainty (see [6] for more de- 
tailed discussions on q 2 dependence). This problem of 
the unphysical q 2 dependence is recently solved by the 
new technology, the so-called QCD factorisation (QCDF) 
approach [13]. Our result is consistent to the one using 
QCDF within the errors caused by the coefficients with 
the higher Gegenbauer terms [14]. The small imaginary 
parts in Eqs. (15) and (16) are the strong phases which 
come from the QCD correction terms in the effective Wil- 
son coefficient in Ref. [10]. The different sign for the 
contributions from the Oi and Oi in B — ► rj'Ks (see Eq. 
(9)) gives the minus sign for the coefficients of the RL and 
RR mass insertions in Eq. (16). It is important to men- 
tion that this sign flip is not due to our using the naive 
factorisation approximation and would not be influenced 
by any other QCD corrections. Note that the absolute 
value of the mass insertion (8j iB )2 3 , which is relevant to 
the b — > s transition is constrained by the experimental 
results for the branching ratio of the B — > X$j decay: 
\(4l,rr)23\ < 1, \(5 d LR , RL h 3 \ < 1.6 x 10- 2 [15]. 

We found that the coefficients for the each mass in- 
sertions are almost the same in B — > (f>K and B — > 
rjK apart from the signs. Accordingly, we shall re- 
parametrise these ratios as: 



J^SM 



4>K S 



-4 s 



/Ism / 



=fye Wl2 e** = 5 L e iBZ * SL + 5 R e iaisSR , (17) 
sifye^V 9 *' ~ 6 L e iaTs5L - 8 R e iaTs5R (18) 



where Om^) and 5\ 2 are CP violating and conserving 
phase differences between SM and SUSY, respectively. 
Sl and 5 R include the contributions proportional to the 
mass insertions (5^ L lr)^3 an d (d RR R l) 23 , respectively. 
Using these parameters, the mixing CP asymmetry is 
given as [6] 



~><f>K s 



S V 'K S — 



sin 2(3+2R 4> cos S 12 sin(6» + 2/3) + R\ sin(2(9 + 2/3) 

1 + 2i?0 cos 812 cos 9$ + R 2 
sin 2/3+2R n > cos 8 12 sm(6 v > +2/3) + R 2 , sin(26»^+2/3) 

1 + 2R ri > cos 8\2 cos 8 n i + R 2 , 



where we use sin 2/3 = 0.73 in our analysis. As can 
be seen from the above formulae, the strong phase en- 
ters only as cos £12 and the small strong phases found in 
Eqs. (15) and (16) lead to cos(5i 2 = 0.99. Thus, we use 
cos 812 = 1 in the following. 

Here let us recall our main conclusions on S^Ks m 
Ref. [6]. function of 9^ behaves as a sin 9^ curve 

taking the value S^Ks = 0.73 at the origin and bounded 
above by 1. A typical behaviour of S^Ks with — 0.5 
and cos<5i2 = 1 is shown as the solid line in Fig.l. In 
the following, we will use this result as a reference and 
fix R<f, = 0.5 and also focus on the region — 37r/4 < 9^ < 
—n/2 where S<f,K s becomes negative. 

Now let us discuss the B — > rj' Ks process and see if 
we can explain out the puzzle of the observed mixing CP 
asymmetries: S^Ks < while S^Ks 1z,Sj/iI>k s - First, 
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-n -n/2 7t/2 n 



FIG. 1. St/,K S (solid line) a,n<fs v i Ks versus 9^, for i?^, = 0.5 
(dashed line for Case 1, dotted line for Case 2, dash-dotted 
line for Case 3 and dash-double-dotted for Case 4). 

we shall show our result without including the contribu- 
tions from these new mechanisms suggested to enhance 
the branching ratio of B — * r/'K in order to see explic- 
itly the different behaviours of S^k s an d S^Ks due to 
the minus sign in Eq.(18). Having some possible SUSY 
models in our mind, we perform a case-by-case study in 
the following. 

Case 1) \6 R \ » \S L \ 
Eqs. (17) and (18) lead to 

R+e** = \5 R \e lllrsSR , (19) 
i^, e i( V = |£ fl |e l(arg,5 * +7r) . (20) 

The CP asymmetry S v 'k s as a function of aig8n{= 9^) 
is shown as a dashed line in Fig.l. \5r\ is fixed to have 
R(j> — \Sr\ — 0.5. As can be seen from this figure, S v 'k s 
is always larger than the experimental data in Eq. (3) 
where S^Ks is within the experimental range. Note that 
the \Sl\ dominated models give a same curve as S^Ks- 

Case 2) \S L \ = |* K | 

In this case, Eqs. (17) and (18) are reduced to: 

\n 

R^K^ 6 * = 2\S L \cos—e l(ais5L+ais5R)/2 (21) 
R n , Ks e w V = 21^1 sin — e ^sS L +^ g s R +^/2 ^ 

where A9 = argjj, — arg^R. We depict S v 'k s as a func- 
tion of (arg^L + arg<5fl)/2(= 6$) for A9 = tt/10 as the 
dotted line in Fig.l. We fix \5l\ so as to have R$ = 0.5. 
The 7r/2 shift appearing in Eq. (22) can be clearly seen in 
the plot. It is also remarkable that in this case, not only 
the phase shift between 9^ and 9^ but also the amplitude 
difference which is given in terms of A9 differentiate the 
behaviour of S^k s and S^Ks ■ I n particular, for small 
A9, no matter what the value of \8l\ is, S^Ks takes a 
value close to sin 2/3. 

Case 3) argdi = argj^ 
In this case, we have 

R*e i6 * = (\S L \ + \6 R \)e iar ^ (23) 
Rr, f e ie n' = A\5\e iaTgSL (24) 

where A | (5 1 = \6l\ — \Sr\. We show our results for S v 'k 3 in 
terms of arg(^(= 9^) in Fig.l for R^ = \5l\ + \Sr\ = 0.5 



and A|5| = 0.2 (dash-dotted line). We found that the 
experimental bound gives a constraint of < A\S\ < 0.4. 

Case 4) arg^R = arg(5i + tt/2 
In this case, we have 

fye"* = V\S L \ 2 + \6 R \ 2 e l ^ sS ^ (25) 
Rr,<e^' = vW+We'^^"' (26) 

where tana = |5r|/|<$l|- In Fig.l, we plot the result 
of S v >k s as a function of arg<5^ + ot(= 6$) for R^ — 
\/\Sl\ 2 + \Sr\ 2 = 0.5 with a = 5tt/4 (dash-double-dotted 
line). With the phase shift of 2a, one can have both S^Ks 
and S n 'K s within their experimental range. 

We should comment that the above model indepen- 
dent analysis can be realised in well known SUSY mod- 
els. For example, the SUSY models with Hermitian fla- 
vor structure which provide an interesting solution for 
the SUSY CP problem [16], have 6 LR = (S RL )* with 
negligible (5ll,rr), which is a realisation of Case 2 with 
A9 = ir. Also in the SUSY seesaw models which are mo- 
tivated by neutrino masses, Srr is much larger than Sll 
[17]. Therefore, Case 1 can accommodate these models. 

As mentioned, another large discrepancy is observed 
in the branching ratio of the B — > rf K process [9]: 

Br e xp . (B ^ ^ _ (55 +i9 ± 8 ) x 1Q -6 ( 27) 

which is 2 to 5 times larger than the standard model cal- 
culation [18]. Since such a large deviation is observed 
only in the B — > rj'K process, the new mechanisms 
based on the peculiarity of rj meson, for instance intrinsic 
charm [19] or gluonium contents of rf [20], have been in- 
vestigated. We shall discuss in the following: the SUSY 
effects to the branching ratio and the impacts of those 
new mechanisms on the mixing CP asymmetry S„/k s ■ 

Let us first discuss the SUSY contributions and also the 
uncertainties from various SM parameters. In general, 
the SUSY contributions can be written as 

Br(B -> r)'K) = Br SM (B -► r)'K) x 

[l + 2cos{9 v , -S 12 )R V , +R 2 ,] (28) 

Note that this equation can be applied to B — > <j>K by 
replacing the indices. The input parameters which are 
used in our above analysis lead to Btsm(B — » rj'K) — 
13 x 10~ 6 . In fact, this value is sensitive, especially to 
the s quark mass and the value of q 2 in our calcula- 
tion. For instance, m s = 0.08 GeV and q 2 = ml/2 
give Btsm(B — > T]'K) = 36 x 10~ 6 . However, such a 
small value of m s enhances the branching ratio of some 
similar processes like B — > irK [18] and also a larger 
q 2 is disfavoured by analysis of k s [6]- A maximum 
enhancement from SUSY contributions can be obtained 
by 9,y — nn, n — 0, 1... and cos<5i2 = 1, which lead to 
Br(B -► rj'K) = 2.25 x Br SM (B -► rj'K) for R n > ~ 0.5. 
Interestingly our solution to reproduce the experimental 
result of S<j,Ks an d S v >k s requires a shift between 9^ and 
9rqf , which may suppress the leading SUSY contribution 
to the branching ratio for B — > <f>K, Thus, it is possible 
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to enhance B — > j/Jf without changing the prediction 
for B — > too much. On the other hand, the other 
similar processes such as B — > 7rA require more atten- 
tion. Apart from its tree contributions, B — » 7rA obtains 
as large SUSY contributions as £? — * r/A. Therefore, we 
must not ignore the limitation given by these similar pro- 
cesses, which would be revealed as soon as more precise 
experimental data from those processes will be available. 

Now we turn to the new mechanisms proposed to en- 
hance Br(B — > rj'K) and its impacts on S v 'k s - We 
rewrite the amplitude in the following way, 

A{tfK) = A^f Ks + A*?g g + G SM + G SUSY (29) 

where G SM and G SUSY are the new mechanism contri- 
butions to SM and SUSY, respectively. Flowingly, the 
branching ratio including the contributions from both 
SUSY and new mechanisms is modified to 

Br(B -» n'K) = Br SM (B -» r/K)(l + r) 2 x 

[1 + 2 cos(0;, - S 12 )R' V , + i?; 2 ,] . (30) 

where r = G SM /A^ s and R' n ,e<' = + 
G SUSY )/(yl^ s +G SM ). Note that Br SM {B -» 77' A) does 
not include the new mechanism contributions. Having 
the gluonium contributions in mind, we parametrise the 
SUSY contributions from new mechanism as: 

~QS^T = a + (^fl)23] + 6 [(*f*) 23 + (4l) 23 ] . 

where (^£j,(x,m)23 and (<^r(rj,))23 have a same coeffi- 
cient due to the penguin process and also a same sign 
since the amplitude is proportional to only the B — A 
transition form factor. Thus, Eq. (16) is modified to: 

,V<-(^)(^ + (^)(^, 

Although the quantitative estimation of r is difficult at 
the moment, the parameters a and b could be computed 
for a given mechanism. For the intrinsic charm con- 
tribution, we have a = b — since it come from a 
tree diagram. For the spectator gluonium contribution 
(G SUSY /G SM ~ V2/(V tb V t *G F ) G| DSY /G| M ), we obtain 
a = —1.2 and b — —585 at the LL order. The spectator 
gluonium process means that the weak b — > sg transi- 
tion (chromo-magnetic operator O g ) accompanied by one 
gluon emission from spectator is followed by two gluon 
fusion into gluonium in rj [21,22]. Using these values in 
Eq. (31), we find that as r increases \6l\ v ' is reduced 
and |<5r|,,< is enlarged. In fact, this does not disturb our 
previous explanation for the discrepancy between S^Ks 
and S V >K S especially because the sign in front of |<5r|<t5 
and |<5fl|^' remain different, which was a crucial point. In 
Fig. 2, we show the result for the branching ratio versus 
Sr)'K s f° r Case 1 — 4 including the spectator gluonium 
contribution. We fix R^ = 0.5 and 6$ = — 5tt/8 in order 
to have S^Ks — —0.2. As can be seen from this figure, 
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FIG. 2. A case-by-case study for the branching ratio of 
B — > jy'A versus the mixing CP asymmetry S v ik s - We 
assume that ~ 0.5 and ^ ~ — 57r/8 which lead to 
S^Ks — —0.2. The parameter r represents the spectator glu- 
onium contribution in SM. 

we can have both of the CP asymmetry of B — ► r/Ag 
and its branching ratio within the experimental limits in 
a significant range of SUSY parameters space. 

To conclude, we have considered possible supersym- 
metric contributions to the CP asymmetry S^Ks an( ^ 
Sr/'Ks ■ We showed that the discrepancy between their 
measurements can be naturally resolved by considering 
the different parity sensitivity of these processes to the 
SUSY contributions from R-sector that lead to R v > < R<p 
and/or a phase shift between 8 V ' and 6^. We also stud- 
ied the observed large branching ratio of B — ► rj'K. We 
have considered the new mechanisms proposed to en- 
hance Br[B — > rj K) and their impacts on S<f,K s ~ SrfKs 
correlation. We have shown that a simultaneous so- 
lution for discrepancy between the CP asymmetries of 
B — > 4>Ks and B — > 77' As and the puzzle of the large 
branching ratio is possible for some SUSY models. More 
prices experimental data would allow us to draw more 
definite conclusions and shed light on compatible SUSY 
models with this solution. 
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